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Abstract 

This is a short review of methods and results of calculations of fission barriers and fission half- 
lives of even-even superheavy nuclei. An approvable agreement of the following approaches 
is shown and discussed; The macroscopic-microscopic approach based on the stratagem of the 
shell correction to the liquid drop model and a vantage point of microscopic energy density func¬ 
tionals of Skyrme and Gogny type selfconsistently calculated within Hartree-Fock-Bogoliubov 
method. Mass parameters are calculated in the Hartree-Fock-Bogoliubov cranking approxima¬ 
tion. A short part of the paper is devoted to the nuclear hssion dynamics. We also discuss the 
predictive power of Skyrme functionals applied to key properties of the hssion path of It 

applies the standard techniques of error estimates in the framework of analysis. 

Keywords: Fission barriers; hssion half-lives; selfconsistent methods; 
macroscopic-microscopic method; Skyrme functional; Gogny functional 


1. Introduction 


Considerable progress in the experimental s\mthesis of heaviest nuclei has been achieved 
recently by the Flerov Laboratory in Dubna and was partially con- 

hrmed in the laboratories at GSI |[il[il[Ii[3[3l, LBNL lUTnand RIKEN lUl. Nonetheless, 
the fundamental question of what is the largest possible atomic number of an atomic nucleus 
is still unanswered. For example, at the time of writing this manuscript the heaviest identihed 
super-heavy element is Z = 118, but even larger elements are not yet excluded. Due to substan¬ 
tial Coulomb repulsion for super-heavy nuclei the stability decreases rapidly when going up in 
system size. Only shell effects allow super-heavy (SH) nuclei with Z > 100 to survive. There are 
two dominant processes of disintegration for such heavy systems: alpha particle emission and 
hssion. It is the subject of this contribution to provide an overview of the state of the art relating 
to the latter process, hssion of SH nuclei. Two different microscopic methods for estimating 
the stability of SH nuclei are used. There are, on the one hand, the self-consistent approaches 
using effective interactions as, e.g., the Hartree-Fock (HF) plus BCS (HFBCS) or Hartree-Fock- 
Bogoliubov (HFB) method, and, on the other hand the more phenomenological macroscopic- 
microscopic (MM) method. The latter is less demanding and thus allows even extensive studies 
of multidimensional hssion landscapes. Although these models differ quantitatively, they agree 
in predicted topological properties as, e.g., the prolate deformed SH nuclei with Z = 100 - 112 
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and N< 170, which are confirmed experimentally for nuclei around ^^"^No in and spherical or 
oblate deformed systems with N= 174 - 184 (see m 21]). 


, Lublin- 
- energy 


The height of the fission barrier, Bf, is one of the most important ingredients to calculate the 
survival probability of SH elements synthesised in heavy-ion reactions. It allows to estimate the 
competition between the fission process and particle emission. To evaluate B /, one needs the 
energy as function of collective deformation, the potential energy surfaces (PES). Ideally, one 
has to figure out all competing minima and saddle points on multidimensional maps. Having 
preconceived ideas about the expected fission landscape, one usually one restricts the search to 
a one-dimensional fission path. In any case, the PES is the basic ingredient for the evaluation of 
fission half-life (HE). 

The last review of theoretical properties of SH elements covering MM and self-consistent 
models was published eight years ago 12211 . Since then substantial progress has been made in the 
development of models and methods es pecially by Lublin-Madrid collaboration 123 2‘ 

Oak Ridge-Warsaw group j 3 11, 23, 28, and Erlangen-Darmstadt team M 

density functional (EDE) HEB approaches - and in the MM model by the Warsaw 132LI33LI34|] . 
Los Alamos school IHl and the relativistic mean-field (RME) IHl . 

In this paper, we are discussing methods to obtain the PES, barrier heights, mass parameters 
and fission HLs using modern approaches like the shell correction method in a multidimensional 
space of deformations and/or the method of energy density functionals (EDE) of Skyrme and 
Gogny types treated self-consistently in the HEB framework. The main purpose is a comparison 
of results of these different approaches and the discussion of uncertainties that affect the data in 
the process of calculation. 

The paper is outlined as follows: Section 12.11 deals with MM methods and data. In Sec- 
tion l2.2l we deliver the discussion of microscopic Skyrme’s and Gogny’s functionals approaches 
and results. The short Section l273l tackles the issue of mass parameters and their effect on the 
HLs. Section [prelates the calculated data to the experimentally known cases. Section [3T2l deals 
with an influence of dynamical treatment of fission process based on the action minimisation on 
fission half-lives. Section|4]discusses an estimate of extrapolation errors for the various observ¬ 
ables. Conclusions are given in Section|5] 


2. Predictions for barriers and life-times 

2.7. Macroscopic-microscopic (MM) approach 

The MM model composes the total energy of the nucleus from two complementary parts. One 
is the macroscopic energy, usually some version of a liquid drop (ED) energy and the other the 
shell correction which describes the influence of single particles (SP) energies and their quantum 
shell structure beyond the ED. Both parts are modelled independently. The ED part is taken from 
the extensive fits within the ED model. The results discussed here were obtained within a MM 
model based on the single particle (SP) spectrum of a deformed Woods-Saxon potential 1^ . The 
SP spectrum obtained from the deformed Woods-Saxon potential is used as input to compute the 
Strutinsky shell correction energy Eor the macroscopic part we used the Yukawa plus 

exponential model 0. The pairing parameters and three parameters of the macroscopic energy 
formula were determined in 14111 by a fit to masses of even-even nuclei with Z > 84 and A > 126 
as given in Id^l . 

In Pigure[T]we display the height of the first fission barriers Bf for these nuclei calculated 
within the MM model as the difference between the lowest saddle point energy and the ground 
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state one Ej. It can be seen that in the whole region of considered nuclei the barriers are smaller 
than 7 MeV. The highest values are obtained for the nuclei ™108i62, and around the 

nucleus Z « 100, =» 150. 
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Figure 1: Contour map of calculated fission baiTier heights Bf for even-even superheavy nuclei. 


Another important observation concerns the behaviour of fission barrier heights with increas¬ 
ing proton number of these heaviest nuclei. One can recognise that the quite high barrier for 
^^®118 i 78 rapidly decreases reaching a value 1.43 MeV for the ^^^126i86 nucleus. Obviously 
such a compound system does not have any chance to survive against hssion in the MM model. 
This is in a contradiction to self consistent models, according to fission barrier for neighbouring 
nucleus is about 12 MeV 143.441. 

Fission barriers for even-even SH elements lf33l obtained within the Woods-Saxon (WS) 
model have characteristic maxima at Z = 108, N - 162 (deformed magic shells) and Z — 114, 
N - 178 (not 184). For the latter system the barrier reaches 6.34 MeV and then decreases with 
Z. This is in contrast with Skyrme SkM* predictions lE^ of the highest barriers for Z = 120 
1 2511 . which is related to the proton magic gap. One has to emphasise that the barriers from the 
WS model are in good agreement with data for even-even actinides; for 18 first and 22 second 
barriers the rms deviation equals 0.5 and 0.7 MeV, respectively ll^l4^ . 

Potential energy surfaces (PES) are often calculated and projected on the triaxial (J5 cos y, 
P sin y) plane. Examples of energy maps in {j3 cos y, j3 sin y) plane, necessary to evaluate fission 
barriers are shown in Eig. |2(The total energy is normalised in such a way that its macroscopic 
part is equal zero at the spherical shape of a nucleus.) 

Calculations of ERDM presented in 14^ confirm aforementioned findings of the WS model 
e.g. region of the highest fission barrier around ^^°Hs and ^^^Lvng, decrease of stability between 
these two maxima as well as beyond N - 184. Nevertheless some qualitative discrepancies 
between two MM approaches can be found. In ERDM barriers are systematically higher by over 
2 MeV. Eission barriers of odd isotopes are also calculated and discussed in the cited paper. 
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Figure 2: Energy surfaces: E - E^^cUphere) for the three of experimentally detected evaporation residues and one 
hypothetical heavier system. 


2.2. Self-consistent theories 

The MM methods with their two independent ingredients (LD model for bulk parameters and 
WS model for SP properties) lead to a rather large number of theoretically independent param¬ 
eters separately fitted in LD and SP parts of the energy. Contrary to this empirical treatment, 
self-consistent approaches start from one density functional model for the energy of a nucleus as 
a whole depending on a couple of densities (local density, kinetic energy density, spin-orbit den¬ 
sity, ...). These total energies can be generated in the HF or HFB framework from given effective 
interactions or modeled directly as density functionals. Such approaches exist in both regimes, 
non-relativistic as well as relativistic. In the following we describe the FIFE approach based 
on two most commonly used (non-relativstic) functionals: the Skyrme type with zero-range nu¬ 
clear forces and the Gogny type where the forces are of mixed nature and both zero-range and 
finite-range components are taken into account. 


4 

































N 

0 

JZi 

E 

c 

c 

o 

^—1 

o 

D. 



148 152 156 160 164 168 172 176 180 184 188 



10 

9 

8 

7 

6 

5 

4 

3 

2 


Neutron number N 


Figure 3: First fission banier (MeV) of SH nuclei in HFB SkM* model. Cross-hatched squares represent observed nuclei. 


2.2.1. Skynne type functional 

The energy in Skyrme-Hartree-Fock (SHF) or Skyrme-Hartre-Fock-Bogoliubov (SHFB) the¬ 
ories of nuclei reads 

£ = (1) 


where p is the local density, t the kinetic energy density, J the spin-orbit density, and so on. (see 

e.g., 00 .) 

There is a variety SHF energy density functionals which differ mainly by their parametrisa- 
tions. One which is often used in the description of the fission process is SkM* 10 . Heights of 
the fission barriers known from experimental data and spurious rotational energies in the ground 
states of even-even nuclei were taken into account in the fitting procedure. Thus, fission proper¬ 
ties should be correctly calculated. 

The rather involved calculations of barriers and the fission HLs of SH elements on a global 
scale were possible and done only since very recently (see e.g., lilSIll). The goal of the 
present section is to recalculate the results of previous estimates for HLs taking into account 
two facts: First, as it was mentioned in many places (see e.g., El 10] where the rotational 
moment of inertia - a part of inertia tensor - is discussed), the Inglis cranking mass is too small 
in comparison to the full ATDHFB one. In this paper, we use the cranking mass but scaled with a 
factor 1.3 to simulate the ATDHFB mass. We follow this principle in the present section. Second, 
a key ingredient is the energy Eq of the collective ground state; in earlier papers we used rather 
crude estimates for Eq and in the following we use an improved semi-classical recipe (Bohr- 
Sommerfeld quantisation) consistent with WKB approximation which is used to calculate the 
tunnelling probability for spontaneous fission. Note that other calculations proceed diifferently in 
those details. For example, reference 0311] used fully quantum mechanically computed Eq withint 
the Bohr Hamiltonian approach and the full ATDHF masses. Figure[2shows the systematics of 
the bander heights calculated with SkM* in the region of Z = 108(2)126 and N - 148(2)188. 
In the experimentally known region, the banders are of the order of 6-7 MeV. The baniers of 
heaviest nuclei are relatively high and reach 8-10 MeV. From the point of view of the bander 
height the most stable nuclei with respect to fission decay are close to Z = 120 - 122 and 
N = 178 - 180. There is a region of low baniers height near the line A = 282. This would 
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predict a valley of fission instability in this particular region of superheavies. On the other hand 
there is experimental evidence ifs^l which contradicts calculations. 
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Figure 4: Ground state energies from WKB quantisation rule (see Eq.[3. Different deformation regions: prolate, y-soft 
and spherical are approximately depicted. Cross-hatched squares represent observed nuclei. 


The calculation of main ingredients of the theory, namely the PES and the collective mass 
parameter were described elsewhere 3. We map the deformation landscape by adding a 
quadrupole constraint stretching the system successively to scission. The PES is then obtained 
as a function of the mass quadrupole moment defined as the average value of the quadrupole 
operator q 2 — {Q 20 ) - (ZiiLi ~ ~ yf))- Two further constraints, namely the triaxiality and 

reflection asymmetry, were applied in the region of the first barrier and beyond. A detailed dis¬ 
cussion is given in |25, 2^. Here we only discuss some new aspects concerning the spontaneous 
fission process. 

The ground state energy measured relatively to the minimum of the potential energy is crucial 
in determining the HLs of the nucleus with respect to fission. The approximate fission probability 
P stemming from WKB theory involves mass parameters B, potential energy V and the energy of 
the collective ground state of the system Eq. The requirement that the phases of WKB functions 
coincide (modulo utt) in the region of their overlap leads to the quantization condition. Then the 
following relation between the quantities involved should be fulfilled for a stationary state 


Sah = (n + l/2)7r. 


( 2 ) 


where, in general, the action Sab (in units of h) calculated between two points a and b is 

r* 

SabiEo)^ ^J2B{q)[Eo-V(q)]dq. (3) 

Ja 

The coordinate q is the leading collective coordinate towards fission {e.g., q = {Q 20 ) q - p), 
n - 0,1,2 etc. In Eq. |2]a and b are the classical turning points determined from the equa¬ 
tion V(q) - Eq solved in the vicinity of the ground state. The Eq. (O is similar to the Bohr- 
Sommerfeld quantisation rule (see e.g., ifsilB . It can be easily solved for each considered nucleus 
with respect to Eq. While the fission probability P strongly depends on Eq the results of calcula¬ 
tion will be different from those where one assumes a constant Eq - 0.5 MeV for all nuclei (see 
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e.g. |56, 5j,]) or quantum mechanically exact values m. The values of the ground state energies 
Eq calculated from the Eq.|2]are shown in Figure 01 It is interesting to observe correllations of 
this quantity with the corresponding ground state deformations . 

For the mass parameter we aim at the adiabatic time dependent HFB (ATDHFB) cranking 
approximation. As it was already mentioned, the Inglis cranking mass parameters are on the 
average smaller then full ATDHFB ones. Therefore to do more realistic calculations of fission 
HLs we assumed according to other authors the scaled fission mass tensor component. The 
scaling factor 1.3 was applied. This leads in general to larger values of Eq [see Eq. (O] but HL 
remains the same as mass parameters are larger. 

The next Figure |5] shows the fission HLs calculated in this way with SkM* parameter set. 
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Figure 5: Fission HLs (sec) of superheavy nuclei in microscopic HFB SkM* model. Cross-hatched squares represent 
obsei'ved nuclei. 
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Figure 6: Total HLs (fission + a) of superheavy nuclei in HFB SkM* model. Cross-hatched squares represent observed 
nuclei. 


The WKB quantisation rule which delivers the ground state energy Eq for fission works well at 
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least in the region of known elements as e.g., in the case of fermium isotopes 115811 . Its extrapola¬ 
tion to SH elements seems to be better than an ad hoc procedure which was assumed previously 
where the ground state energy Eq was equal to 0.7 where ^zpg is zero point energy, or it 

is a constant equal to 0.5 MeV. 

After taking the a-decay into considerations the HLs for some nuclei change substantially. 
To estimate or-decay half-lives, we used the standard Viola-Seaborg expression 15911 with the pa¬ 
rameters from Ref. I^oll . The resulting HLs are shown in Figure|6] According to our calculations 
the longest total HLs correspond to nuclei in the vicinity of Z = 112 and N - 182, 184. The 
region of most stable nuclei is located in the vicinity of Z ~ 112, A ~ 182, 184 ll^ . It is worth 
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Figure 7: Regions of symmetrically and asymmetrically fissioning nuclei in HFB SkM* model. The filled circles (blue) 
represent nuclei decaying mainly by alpha emission. The isolines show approximate total HLs (in seconds). Experimen¬ 
tally registered nuclei are depicted as thick (green) squares. 



to mention the existence of two main modes of hssion. The symmetric mode (the same spherical 
fragments) which dominates in the region of light SH nuclei and the asymmetric one prevail in 
the region of still heavier systems. Both symmetry regions are shown in the Figure|2l The sep¬ 
aration line is close to A ~ 280. The asymmetric hssion region is concentrated at Z ~ 118 and 
N ~ 178. It is observed that the heaviest nuclei {N > 186) decays probably in the asymmetric 
way. The thick squares shown in Figure [T] correspond to nuclei for which experimental data for 
hssion and/or alpha decay HLs are known ifssll . These trends are in agreement with the survey 
of if^ using a couple of other Skyrme parametrisations. 

Fission barrier heights and HLs with respect to spontaneous hssion as obtained from self- 
consistent HFB SkM* calculations differ considerably from the ones presented in the previous 
section where MM models were used. The main difference are the barriers. They are usually 
higher in the case of HFB SkM* model up to 2-3 MeV at some cases. This feature is common 
for all non-relativistic self-consistent approaches. 

There are of course discrepancies between the theoretical and the experimental data, the 
largest ones in the region of the middle mass SHE. Since those systems belong to the region of 
shape coexistence some further increase of spontaneous hssion HLs is anticipated due to shape 
mixing. In order to improve upon the theory of hssion we shall consider the inclusion of several 
collective coordinates and the dynamics of spontaneous hssion process, the improved energy 
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density functionals 116ill , and the full ATDHFB inertia parameters. The work is in progress. 
Some introductory results have been already published 0291] . 


2.2.2. The Gogny functional 

The Gogny force ll62[l is a non-relativistic effective phenomenological interaction widely used 
to describe the low energy dynamics of the atomic nucleus. It is intended to be applicable to all 
possible nuclei from drip line to drip line using a limited set of parameters adjusted to global 
properties. In addition to the standard central, spin-orbit and Coulomb potentials it contains a 
phenomenological density dependent term depending on the density raised to the power y = 1 /3 
introduced to incorporate the saturation mechanism of nuclear forces. In many aspects, it is sim¬ 
ilar to the Skyrme 16311 interaction sharing with it the spin-orbit potential. Coulomb and the idea 
of a density dependent repulsive part. However, Skyrme’s zero range and gradient terms present 
in the central potential are replaced in the Gogny force by a hnite range potential linear combi¬ 
nation of two gaussians. The main advantage of the hnite range is that the corresponding ma¬ 
trix elements are free from ultraviolet divergences allowing the treatment of the particle-particle 
(pairing) channel with the same interaction. 

In most of the hssion studies carried out with the Gogny force, the focus has been in the 
characterisation of just one of the hssion observables, namely the spontaneous hssion lifetimes 
(or derived quantities like hssion barrier heights and widths) in even-even nuclei. Studies aimed 
at the description of induced hssion in the framework of hnite temperature HFB are scarce ll^ 
and deserve more attention. This is the same as with hssion in odd mass nuclei where the un¬ 
derstanding of the specialisation energy is of fundamental interest. In this respect, there is only 
one example with the Gogny force ll6.5ll . Recent advances in the experimental techniques allow 
for the study of the evolution of hssion parameters with spin ll66ll . The expected decrease of the 
barrier height with spin and the evolution of triaxiality in the hrst barrier has been studied with 
the Gogny force in Ref 0. 

At this point, it has to be mentioned that the ability of any of the Gogny parametrisations 
to reproduce hssion observables (mainly spontaneous hssion HLs) is spoiled by the strong de¬ 
pendence of those quantities with the value of the collective inertia in the forbidden region. As 
the collective inertia strongly depends on the pairing gap as the inverse of its square it is to be 
expected that pairing properties should play a central role in the theoretical description of hs¬ 
sion. In this respect, we can mention two studies ^ with DIS, DIN and DIM where the 
spontaneous hssion lifetimes of the neutron rich uranium and plutonium isotopes have been con¬ 
sidered. The variation of the HLs with the pairing strength was found to be enormous in absolute 
value although the trend with neutron number was roughly preserved. As a consequence of this 
variability, the results also depend upon the principles involved in the determination of the path 
to hssion. The recent results of ll70ll using the minimum action principle with pairing degrees 
of freedom are markedly different from those using the minimum energy principle {see Section 
13.2b . The reason again is the strong dependence of the collective inertia on the pairing gap (it 
goes as the inverse of the square of the pairing gap) that has a strong impact on the action. 

The PES calculated in a microscopic model like HFB theory provide a lot of information on 
the hssion properties of SH elements. HLs, preferred decay channels and hssion types (sym¬ 
metric, asymmetric, highly asymmetric) can be deduced by the analysis of PESes. As it was 
mentioned above, results derived within HEB theory are obtained with some systematic theoret¬ 
ical errors {see Section |4]l. Nevertheless, it is reasonable to compare results between particular 
isotopes and distinguish groups of nuclides with specihc decay properties. An early stud y of SH 
region in the HEB model has been performed by Bruyeres-le-Chatel group in 171.72. 73 .HItI 
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Figure 8: Axial fission barriers of selected isotopes across the SH nuclei. Note that in the original paper (U the 
quadrupole moment Q 2 was a factor of two smaller. 


and Lublin-Madrid group in S. The most thorough investigation using the DIS Gogny in¬ 
teraction in the axial symmetry regime has been published in 11241] . Details of the method are 
described in 115ill . 

Fission characteristics are determined not only by the height of the fission barrier but also 
by its shape and localisation in the deformation space. From the general overview of the fission 
barrier shapes across the SH nuclei one can find that the number of protons affects mainly the 
height of the barrier whereas substantial changes of the barrier profile take place along isotopic 
chains. In consequence, fission properties depend mostly on neutron number. Selected fission 
barriers of SH isotopes calculated in the axial regime are presented in Fig. [8] 

The lightest SH elements, with neutron number below N - 170, are well prolate deformed 
in the ground state. The saddle point in these nuclei has very high energy, often exceeding 8 
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Figure 9: Logarithm of the spontaneous fission HLs in seconds computed with the HFB model and the DIS Gogny force. 
Nuclides with dominant fission are represented by filled symbols whereas those with dominant n decay are marked by 
open symbols. Dotted line indicates isotopes which fission barrier is presented in Fig. [8] 
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MeV. Going along any isotopic chain in this region from N = 150 to = 162 one can observe 
increasing width of the fission barrier while the barrier height remains almost unchanged. As a 
consequence, the action integral across the barrier increases and longer HLs are observed (see 
Fig. |9]). The centre of the region of extended stability is predicted for ^^°Hs with the longest 
fission HL log(fsF/s) = 12.83. The corresponding nucleon numbers, namely N - 162 and 
Z = 108, are sometimes called “deformed magic numbers”. The special characteristics of the 
A = 162 and Z = 108 system are correlated with the decrease of the neutron pairing energy in the 
ground state and a small value of the proton one 11241] . This is consistent with the aforementioned 
MM and Skyrme model results. 

In most of the isotopes below N - 162, the fission barrier goes through reflection symmetric 
shapes. Only in the nuclei in the transitional region between the actinides and the SH elements 
with A ~ 250 octupole shapes play an important role in the determination of the fission barrier 
height. The reflection symmetric character of the barriers brings about symmetric fragment mass 
distribution in this region. The barrier height may be reduced by a few MeV in the whole region 
if non-axial deformations are included in the calculations EH. 

Beyond N = 158, we And several different trends in the structure of the PES. First, the top 
of the fission barrier splits into two peaks: the first one is lower and the second one higher. In 
isotopes heavier than A = 162 the height of both saddles decreases and, in consequence, HLs 
of subsequent isotopes become shorter. The second, reflection symmetric saddle may be ad¬ 
ditionally diminished if other degrees of freedom characterising the nuclear shape are allowed. 
Octupole deformation cuts off the second barrier, but up to A = 280 such reflection asymmet¬ 
ric barriers extends to large quadrupole deformations and can not become the dominant fission 
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channel. In heavier nuclides the width of the octupole deformed hump in the fission path is com¬ 
parable to the reflection symmetric one and asymmetric fission can be observed as a dominant 
fission channel. Also triaxial shapes can decrease the height of the second barrier but at the cost 
of increasing the inertia and therefore the contribution to the action remains unchanged 13 • 

In the same range of neutron numbers, the ground state deformation decreases and energy 
difference between prolate and oblate minimum reduces. In isotopes with 178 < N < 182 the 
oblate minimum becomes the ground state. When non-axial deformations are taken into account, 
the local minimum with higher energy becomes a saddle point and there is a path along the y 
plane connecting it with the ground state. 

As a consequence of the aforementioned changes of the PESs, fission HLs decrease with 
the mass of nuclides. Minimal values, not longer than fsp = 10^ s, are reached around neutron 
number N - 170 in each SH element. 

The first peak of the barrier keeps the axial reflection symmetric character in the whole region 
above N - 158. Its height shows a growing trend with both neutron and proton number in 
isotopes between N - 170 and N - 182 which is responsible for increasing fission HLs up to 
N - 182. The isotopes with two neutrons more (N = 184) have got features of magic nuclei: 
well pronounced spherical ground state with small pairing energy and large shell gap in neutron 
single-particle energy spectrum. However, a lower fission barrier induces here larger fission 
probability than in lighter isotopes. There is no strong indication of magic character of any 
proton number in this region of SH elements in the calculations with the Gogny force. 

Beyond A = 184 the second, octupole deformed barrier drops away below the ground state 
energy level and fission barrier becomes substantially narrower. This change dramatically re¬ 
duces fission HLs by several orders of magnitude with adding subsequent pair of neutrons. On 
the other hand, in the nuclei with neutron number greater than A = 184 octupole deformation 
decreases the ground state energy by a varying amount that can reach the 2 MeV. This decreas¬ 
ing of the ground state energy work in the direction of increasing again the height of the fission 
barrier. However, this effect is not strong enough as to counteract the above mentioned reduction 
of the fission HLs. 

Two regions of extended stability against fission can be found in the SH elements, see Lig. 
0 The first is around the deformed magic nucleus ^^°Hs where fission HLs reach the level 
of lO'^ s. The second region covers isotopes with neutron number A = 182. The HLs in 
elements with large proton number reach here much larger values extending far above fsp = 10^° 
s. It is important to notice the lower stability against fission in isotopes around A = 170 which 
is presently in the focus of experimental investigations. In nuclei heavier than A = 186, a 
decrease of potential energy with deformation dictated by LD bulk properties dominates over the 
stabilising shell correction. High fission barrier can thus not be achieved and the fission stability 
drops rapidly. 

The main decay process competing with spontaneous fission in SH elements is a radioactiv¬ 
ity. The predicted a decay HLs are shorter than fission ones in proton rich isotones. Thus, it can 
be noticed that fission is the dominant decay channel in three regions 11241] . Lirst, in the vicinity 
of deformed magic nucleus ^^°Hs symmetric fission can be found. Second, asymmetric fission 
is dominant for nuclei elements lighter then LI (Z = 114) and below magic neutron number 
A = 184. Linally, for heavy isotopes with A > 186 and Z < 118 again symmetric fission occurs 
with very short HLs. 

In those isotopes beyond A = 190, the macroscopic part of the nuclear potential makes fission 
barriers too small to allow for sufficiently long HLs that could be detected experimentally. The 
existence of heavier nuclei is possible only in exotic ground state configurations. Bubble nuclei, 
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spherical symmetric systems, without nuclear matter in the centre, are predicted for charges 
Z > 240 and masses A > 700 zS 22. 24, 75]. Lighter nuclei with Z > 170 and A > 450 can 
exist in semi-bubble configuration where central density is reduced to half of the bulk nuclear 
matter density. Another type of exotic shape that may exist beyond SH nuclei is the one of a 
torus 0. The HFB theory with Gogny force predicts a toroidal structure with energy lower 
than spherical configuration for nuclei with Z > 140. The stability of such exotic shapes is 
questionable and it has to be confirmed in further investigations. 


2.2.3. Relativistic Mean Field 

Fission barriers in SH nuclei has been also investigated in covariant density functional theory 
1 ( 3 ^ . Modern parametrisation DD-ME2 and DD-PCl was used together with modified version 
of traditional NL3*. It was found that both triaxiality and octupole deformation substantially 
affect the second barrier in nuclei with Z > 112. It was found that the first barrier remains axially 
symmetric in the whole region. 

Triaxial deformation decrease energy of the saddle point as much as 3 MeV in some isotopes 
with N < 174. The NL3* and DD-ME2 predict the lowest saddle point energies among the all 
models of nuclear forces. The DD-PCl parametrisation gives slightly larger values. 


2.3. Mass parameters and their impact on half-live times 

The collective mass along the fission path influences the fission HLs with the same impor¬ 
tance as the barrier does. This can be seen from formula Q for the action integral where the mass 
B{q) and the potential energy V{q) enter the action as a product. A relatively recent discussion 
of the mass tensor calculations can be found in j^ . 

Mass parameters are usually calculated within perturbative cranking approximation (or Inglis 
inertia) or adiabatic time dependent HP or HPB model where one differentiates the fields of the 
model numerically instead ( see, e.g., I^ls^lsil 31 MM)- 


The shapes of the mass tensor components are fairly similar in both mentioned cases but 
the values of specific components differ substantially. The simple cranking masses calculated in 
e.g., one dimensional space are usually smaller then ATDHP cranking masses. This influences, 
in turn, the fission HLs. In the EDP of the Gogny type models of nuclear structure where the 
barriers are much higher then ones from e.g., MM or RMP models, these small masses lead to 
smaller HLs, which sometimes coincide with experimental ones. On the other hand, the Skyrme 
EDP give smaller barriers then Gogny ones and the ATDHP cranking masses are numerically 
calculated. They generate HLs much larger or comparable to that of the Gogny model. In the 
MM models the barriers are believed to fit rather well experimental data and one uses the simple 
perturbative cranking approximation for the masses. Unfortunately that leads to too small (in 
comparison with experiment) HL times. The other feature is the collective inertias strongly de¬ 
pend upon the pairing strength what makes the pairing degree of freedom an important ingredient 
(see the section on dynamics of fission (I3.2l i). 

The above discussion shows that there is still a multitude of different approaches to compute 
the collective mass along the fission path. This hinders a direct comparison of the final fission 
HL. It is a challenge for the future projects to remove all of these differences in the calculations 
of the nuclear inertia by adopting the same ATDHP cranking model everywhere. Only then it 
will be possible to judge the models and the validity of the barrier heights as well as the validity 
of final HL. 
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3. Comparison of calculated fission barriers to experimental data 

There are two main methods of calculating fission HLs and barriers, namely the static method 
(minimisation of the potential energy and extraction of the saddle points) and dynamic (minimi¬ 
sation of the action integral). Most of the calculations (and these in this paper) concerning fission 
barriers are of the static type. 

3.1. Static fission barriers 



Q20 (barn) 

Figure 10: Fission barriers of ^^^Hs as functions of the quadrupole moment Q 2 Q calculated in MM model (green, filled 
squares), Skyrme HFB approach (blue, filled circles) and Gogny HF model (red, filled triangles). Dashed lines represent 
the raw axial barrier. Solid lines correspond to the barriers with y. The values of Q 2 Q for MM curves are close to the 
actual, but not completely accurate. 


We start with the comparison of the fission barriers for ^®®Hs as an example of mid-mass SH 
isotope calculated in models discussed in the paper. The barriers derived in the MM, Skyrme 
HFB and Gogny HF models are shown in Fig. [10] as functions of quadrupole moment 220- The 
shapes of all PES are very similar and there is relatively good agreement of barrier heights. The 
differences in the peak height (the saddle point) reach 1 MeV. Rather small differences remain¬ 
ing are hard to explain because the models are very different in their structure. The differences 
shown will weakly impact on calculated fission HLs. Fig. [TT] shows four sets of calculated fis- 
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Figure 11: Fission barriers predicted by various models: black - WS model j^, green - FRLDM 1^ . blue - SkM* 1^ . 
red - RMF with NL3 parametrisation (3^. 


sion barriers Bf for even-even isotopes of Z = 112,114 and 120; from the WS model Il32|] . the 
macroscopic-microscopic Finite Range Liquid Drop Model (FRLDM) BSll . FIFE calculations 
with the SkM* interaction |25] and the RMF model with NL3 parametrisation 1^ . Results of 
rather comprehensive calculations were selected; all were made for dozens of nuclei and included 
important triaxial and mass-asymmetric deformations. One has to emphasise that three of those 
models well reproduce experimental barriers in even-even actinides; the SkM* model overesti¬ 
mates actinide barriers by 1-2 MeV. In spite of this, the differences in predictions in Fig. [TTIare 
conspicuous. In order to quantify this, we refer to the known spontaneous fission (s.f.) HLs; 0.8 
ms in ^^^Ds, 97 ms in ^^"^Ds and 130 ms in ^^®F1. The values predicted some time ago from 
the sub-barrier action minimisation based on the WS model were, respectively, 71 ms, 4 s and 
1.5 s (1-2 orders of magnitude overestimate; only axial deformations were included in isHl.but 
axially symmetric shapes at the barrier in these nuclei were actually confirmed in ll32|] L One can 
observe that the FRLDM barriers for Z = 112,114 are much higher than in the WS model and 
could hardly produce HLs in agreement with experiment. As we have checked, this huge differ¬ 
ence in barriers in the two MM models is not related to shape parametrisation. The spontaneous 
fission HLs for the three nuclei calculated in along static path (no action minimisation) are 
(surprisingly, in view of barriers similar to ours) by many orders of magnitude too short. On the 
other hand, RMF calculations produced barriers much lower than in the WS model which 
could hardly explain the measured evaporation residue cross-sections in fusion. 


15 


































In Table [Tl]^ collect some of the theoretical predictions of hssion barrier heights based on 
the FRLDM SHF - self consistent Hartre-Fock method ifs^l with the SLy6 Skyrme in¬ 


teraction Il83n and the Extended Thomas-Fermi plus Strutinsky Integral model (ETFSI) based 
on Skyrme SkSC4 functional ll^ in relation to experimental ones 1851 and our predictions in 
a MM and SHEB(SkM*) models. Note that the lower limits for the hssion barrier heights are 
evaluated in Ref. ll^ . As we can see, experimental and MM calculated barrier heights are 
in agreement while ERLDM, SHE and SHEB generally overestimate the barrier heights. Sur¬ 
prisingly the modern version of the MM calculation presented in Ref. ll35n (including nonaxiall 
shapes) gives values comparable to axially-symmetrical calculations of SHE and overestimate 
the experimental barriers in a similar way. This tendency is visible, to a lesser extent, also for 
heavier systems. Thus for the isotope ^^^Dsne the ERLDM model gives a value of 9.25 MeV for 
the barrier height while the value obtained in our approach (6.22 MeV) is about 3 MeV lower, 
whereas the experimental data indicate 6.4 MeV. Eor the neighbouring isotope ^^^Lvng the exper¬ 
imentally predicted fission barrier is also 6.4 MeV and we obtain also almost the same value 6.28 
MeV, but estimations based on the ERLDM predict a value of 9.46 MeV, almost 3 MeV larger. 
As a consequence, such high barriers result in cross sections for Z=114,116,118 and 120 which 
overestimate the experimental data by several orders of magnitude ll^ 871. The ETESI(SkSC4) 


model significantly underestimates the barrier heights for and ^*®Cni 74 for heavier sys¬ 

tems the agreement is however much better. One must keep in mind that ETESI(SkSC4) and 
SHE(SLy6) models do not include triaxial nor higher order nonaxial degrees of freedom and the 
observed discrepancies probably can be explained by the absence of these variables at the saddle 
point conhguration. 


3.2. Fission dynamics 

As it was shown years ago and also in recent publications |29 the min¬ 

imisation of the reduced action integral instead of the potential energy leads to larger fission 
probabilities or equivalently to shorter hssion HLs. The latter are generally reduced by one to 
two orders of magnitude in comparison to statically calculated HLs. At the same time, the calcu¬ 
lated hssion barriers increase by 1-2 MeV. While these kinds of calculations are rather involved, 
one usually performs static calculations of HLs integrating the action along the potential energy 
valley leading to the hssion of the nucleus. This is common to all calculations. 


The dynamic treatment of some internal degrees of freedom like e.g., pairing gap {see 119 111 
and the very recen t paper if^ l or particle number huctuation parameter /I 2 in Lipkin-Nogami 
pairing model 29, seems to be more signihcant than a similar treatment of space-type 
deformation coordinates described above. While the latter is already well known to be similar 
for all instances the former depends on the considered nucleus. Eor a few cases, it was shown in 
the above cited papers that a dynamical treatment of the pairing held is important and one should 
include pairing Huctuations as a dynamical variable in the calculation of HLs. The dynamical 
change of e.g., the pairing gap being the consequence of action minimisation corresponds to the 
hssion path on which the collective mass parameter is generally smaller than one on the static 
pathway and at the same time the hssion barrier increases as compared to the static one. The 
hnal HE is reduced again. The reduction factor is around one order of magnitude. Eor reasons of 
space, we do not discuss this type of calculations in detail here. 
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Table 1: Comparison of fission barrier heights with another theoretical evaluations: SHF(SLy6) (S^, SHFB(SkM*) - 
this paper, FRLDM l3^ . ETFSI(SkSC4) with Skyrme SkSC4 force fs^l . WS IS - present paper, and experimental data 
taken from ISl . 


Nucleus 

SHE 

(SLy6) 

SHEB 

(SkM*) 

ERLDM 

ETESI 

(SkSC4) 

WS 

EXP 

2S4112i72 

6.06 

4.31 

7.41 

2.2 

4.29 

5.5 

286112i74 

6.91 

4.74 

8.24 

3.6 

5.01 

5.5 

28^114174 

8.12 

6.24 

9.18 

6.1 

5.53 

6.7 

290114i76 

8.52 

7.13 

9.89 

6.6 

5.83 

6.7 

292114i78 

- 

8.14 

9.98 

7.2 

6.34 

6.7 

^‘^^116176 

9.35 

8.72 

9.26 

6.5 

6.22 

6.4 

294i 1 ^ 
iiOi78 

9.59 

8.84 

9.46 

7.2 

6.28 

6.4 


- 

8.54 

9.10 

7.2 

6.07 

6.4 

'‘'"118176 

- 

9.19 

8.48 

6.6 

5.99 

- 

118178 

- 

9.47 

8.36 

7.0 

6.04 

- 

'‘''118180 

- 

9.28 

8.05 

7.4 

5.72 

- 

'‘'^120176 

- 

9.48 

7.69 

6.2 

5.64 

- 

'‘"120178 

- 

10.05 

7.33 

6.6 

5.50 

- 

3001 OA 

- 

9.91 

7.01 

6.8 

5.05 

- 

3021 AA 

- 

9.82 

6.07 

7.2 

4.66 

- 

3041 A A 
iZUi84 

- 

9.54 

4.86 

6.8 

4.20 

- 


4. Uncertainties on barriers and correlations with other observables 


So far, we have mostly presented the results for the various models separately. A direct 
comparison of barriers for four selected models was done in Figure[TT]and it shows that the pre¬ 
dictions on fission properties can vary widely amongst the models. It is the aim of this section 
to explore the variations of predictions on a broader and more systematic basis. To this end, we 
will exploit statistical analysis built on least-squares techniques for optimising parametrisation. 
We will briefly summarise the least-squares techniques, then discuss the various influences on 
barriers by explicit variation of force parameters, and Anally present systematics of extrapolation 
uncertainties and correlations with system properties from statistical analysis. The basic princi¬ 
ples of error estimates and correlations analysis have been explained in great detail in 11921] . A 
practical application in the context of SHF with inspecting a broad range of observables is found 
in i^ . Following the scope if this paper, we concentrate here on fission properties. As we have 
seen in figure [TT] that the isotopic and isotonic trends of fission properties are about the same for 
the different models, we confine the study to one representative test case, fission properties of 
^®®Hs as a prototype super-heavy element. 


4.1. The fission path o/^®®Hs 

The fission path is a set of mean-field states lO^) comprising a range of relevant quadrupole 
deformations q. In self-consistent mean-field models, it is generated by mean-field calculations 
with quadrupole constraint, see e.g. 17811 . This yields the potential-energy surface (PES) for fis¬ 
sion. Figure [12] shows the PES for our test case It is what we call the “raw PES” which 
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deformation energy curve Hs, SV-min 



deformation P2 


Figure 12: Illustration of the potential-energy surface (PES) along the axially symmetric fission path of ^®^Hs and its 
four key quantities: curvature of the PES at the ground state, fission barrier, mass parameter B at the ground state, and the 
action integral = V®(WES^^'Siub) with £sub = '^max“4.2MeV as appropriate average quantities in the tunnelling 

regime. The PES was computed with quadrupole constrained SHE calculations using the parametrisation SV-min. 


emerges as expectation value of energy for given deformation. The actual fission calculations 
will invoke also the mass parameter along the path and augments the raw PES with quantum 
corrections (zero-point energies) 07811311] . This modifies the PES and with it the barriers. Eor- 
tunately, the modifications are moderate for This allows us to confine the discussion of 

trends and coiTelations to properties of the raw PES. Eor simplicity, we will also confine the 
considerations to the fission path along axially symmetric shapes. It is known that triaxiality 
produces for ^®®Hs typically about 1-2 MeV lower baiTiers (cf. Sec. [3] 1 and Eig. [TOl) . However, 
axial shapes suffice for the purpose of comparing different parametrisation and studying trends 
and correlations. 

Eor a systematic analysis, it is advantageous to reduce the information on the PES to a few 
key quantities. Eigure [12] illustrates the choice of these key observables. The most obvious 
quantity is the height of the fission barrier B f. In the following we consider Bf as a difference 
between the saddle point of PES and the energy of mean-field ground state. (Note that the 
effective barrier seen in experiment is rather the difference between the saddle point energy and 
the collective ground-state energy which lies typically 1-2 MeV above the mean-field ground 
state energy.) The collective ground state is characterised by the curvature d^E(q)ldq^\ and 
the mass parameter (not shown here) at the minimum. Another key piece for determining fission 
HE is the WKB action integral (see Eq. O for tunnelling. This is a very involved quantity as it 
depends on entry and exit point for the tunnelling regime which, in turn, depend sensitively on 
the energy of the collective ground state. We want here a simpler measure concentrating on the 
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mass parameter in the tunnelling regime. As this parameter varies considerably along the path 
(see figure [TO we consider an average taken over the core of the tunnelling regime. To that end, 
we consider the action integral S (see caption) in the form as it enters the tunnelling integral. 
In order to concentrate on the effect of the mass B{q), we confine the integration regime safely 
inside the turning points and take a fixed energy of 4.2 MeV. Of course, the actual value for the 
action integral depends on this reference energy. But we are interested in comparing differences 
which develop when varying the forces. The trends seen in such variations are the robust against 
the detailed choice of reference energy and integration limits. 
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Figure 13: The axially symmetric fission barrier of ^®®Hs for a great variety of published models and parametrisation. 
The basic citations for the models and parametrisation are: Gogny I^ .SkM* l4^ . S kP p3l . SkT6 loSl . SLy6 Wh . SkI3 

BSkl (H, SV-min [l^, UNEDF2 fiolll , NL-Z 2 Q 0 I, DD-PC and DD-ME fToak 

Before proceeding tox^ analysis, we illustrate in figure[T3]the broad span of predictions for 
the most important quantity of the fission path, the fission barrier of To that end, we have 

fetched a couple of published parametrisation from early first generation models to latest devel¬ 
opments. There are three formally different models involved, the Gogny force to the left, a large 
group from SHF in the middle, and three RMF models at the right. Even within one model fam¬ 
ily, SHF or RMF, one sees a large variation of results and it is obvious from this figure that many 
of the parametrisation have to be ruled out for fission calculations. For SHF we can establish a 
trend of Bj with the effective mass m* jm of a parametrisation where low m* jm is related to large 
Bf and vice versa. This argument does not apply to the three RMF parametrisation which all 
have about the same low m* jm. Here it is a strongly varying incompressibility K which causes 
the different Bf. After all, the example of figure [T3] demonstrates that an arbitrary collection of 
parametrisations does not allow to conclude unambiguously on the relations between properties 
of a force and predicted observables (here Bf). Statistical analysis based on fits provides a 
well defined strategy to disentangle the various influences and to estimate uncertainties in pre¬ 
dictions. This is what we will unfold now in this section. Before carrying on, we want to remark 
that statistical analysis can unravel trends, uncertainties, and correlations for a given model. Be- 
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sides that we have to be aware of a systematic error which sneaks in by the choice of a functional 
form for the mean-field model. These can be explored by variation of a model and comparison of 
different models as exemplified in ll^ . We will not pursue the quest for systematic errors 
here. 


4.2. Brief review ofx^ analysis 

The Skyrme-Hartree-Fock approximation is a typical representative of presently used self- 
consistent nuclear m odels. Its structure can be deduced from general arguments of low-momentum 
expansion lll04ll05n while the remain ing model parameters are determined by adjustment to em¬ 
pirical data, for reviews see Strategies for such adjustm ent h ave much developed over 

the decades. Present standard are least-squares (x^) fits sta rted i n | l_08l] and steadily extended to 
include more and more data, for recent examples see, e.g., il0dll09l]~ The least-squares (x^) fits 
and subsequent analysis of extrapolation e rrors as w ell as correlations are a standard technique 
and well documented at several places 1 92, 1 id. Ill ]. The tasks is to optimise a set of model pa¬ 
rameters p = (pi, ...,p^p) such that arepresentative set of observables {(9;, / = l...A^^) is optimally 
reproduced. This is achieved by defining a quality function ;t'^(p) from the sum of squared er¬ 
rors and minimising it. Extrapolation uncertainties and correlations between observables can be 
estimated from the behaviour of;^'^ in the vicinity of the minimum. To that end, one defines prob¬ 
ability distribution of “reasonable” parametrisations as VT(p) oc exp(^-^x^(p))- The (statistical) 
expectation value of an observable A is then A = J d^'’pA(p) W(p). Extrapolation uncertainties 

quantify the fluctuation around the mean value and are computed as AA = '\J(a - A^ . The 

correlation between A and B, also called covariance, becomes cab - ~ B^KAAAB). 

The actual evaluation is usually done in terms of Gaussian integrals because x^ can be well 
approximated by a quadratic form in p, see e.g. 


4.3. Choice of fit observables and force parameters 

Eor the fi t observables entering the x^, we consider here the same set as it was used in the 
survey EM- binding energies, key parameters of the charge form factor (r.m.s. radius, diffrac¬ 
tion radius, surface thickness), odd-even staggering of binding energies, spin-orbit splittings in 
semi-magic nuclei selected to be well described by a pure mean-field model. A straightforward 
fit to these ground state data yields the parametrisation SV-min. It embraces all information how 
the given ground state data determine features of the parametrisation and how much they leave 
undetermined. 

A key feature of a model are the nuclear matter properties (NMP). We will consider here 
the four NMP which characterise the response properties of nuclear matter: incompressibility 
K, isoscalar effective mass m* jm, symmetry energy ctsym, and TRK sum-rule enhancement a-trk 
( for a detailed definition see e.g., llOOll i. These NMP have been used since long as benchmarks 
in various nuclear models and thus they have developed an intuitive meaning. Moreover, there 
is a one to one relation between NMP and the force parameters of the volume part of the SHE 
functional. Thus NMP are equivalent to force parameters. Dedicated variation of NMP is useful 
tool to explore influences within a model. Thus we will also consider parametrisations which are 
fitted as SV-min, but with an additional constraint on NMP. There is, for example, the parametri¬ 
sation SV-bas which was fitted while freezing K, m* jm, asym, and /fTRK at values such that one 


20 





















obtains additionally also a good description of giant resonance excitations OlOOll . The differ¬ 
ence between S V-bas and S V-min is one indicator for the influence the these four NMP. We will 
consider in section l4~6l an even greater variation of parametrisations. 


4.4. Extrapolation uncertainties 



0.2 0.3 0.4 0.5 0.6 0.7 

deformation P2 



0.2 0.3 0.4 0.5 0.6 0.7 

deformation P2 


Figure 14: Raw PES (lower panels), mass parameter (middle panels), and ZPE (upper panels) along the axial and 
reflection symmetric fission path of The left panels show results for SV-min, the right ones for SV-bas. Each case 

is shown together with the extrapolation uncertainties from;^^ analysis. 


Uncertainties for extrapolated observables are the simpler signal from analysis. Thus 
we start with a quick glance at uncertainties for the present test case, fission of Figure 

[T4l shows the ingredients for a fission calculation, potential-energy surface (PES), zero-point 
energy (ZPE) and mass parameter, all for motion along axially symmetric quadrupole defor¬ 
mation expressed as the dimension-less q = [52 (for details of the collective picture behind see, 
e.g., OTSilsdl t. The figure label indicates “raw PES”. This is the energy directly from a quadrupole 
constraint SHE calculation. The collective fission dynamics employs the PES corrected by ZPE 
178, ^ . The pattern for the corrected PES are very similar to those for the raw PES and thus we 
have skipped this extra panel. The mass parameter is obtained by self-consistent cranking, also 
called ATDHE cranking |78, All three quantities are shown together with the extrapolation 
uncertainties from;^'^ analysis (see section lATl l drawn as error bars. Note that the PES are scaled 
to zero energy at the minimum. There is thus no uncertainty at this point and only slowly growing 
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uncertainty around it. Sufficiently far away from the minimum, the etTor bars are of similar size 
along the path. The same holds for the ZPE. More fluctuations of uncertainties are seen for the 
the mass parameter. Some points have very large uncertainties and some points for SV-bas did 
not produce a reasonable result at all (that is why there are missing points in the curves). These 
artifacts are due to round-off problems in the very demanding error analysis combined with the 
demanding computation of mass parameters. These problems arise typically in regions where 
many level crossing exist. 

The uncertainties for the unconstrained fit SV-min are larger than those for SV-bas. This is 
plausible because SV-bas is constrained by four NMP which reduces uncertainties. However, 
the reduction is only about factor two. The constraint on K, m* fm, flsym, *'trk thus has only 
partially impact on fission properties. There remains an equally strong influence from the other 
force parameters. The uncertainty of the PES at the barifer is 1 MeV for SV-min and 0.6 MeV 
for S V-bas. This shows the benefit of adding more information from NMP, alias giant resonances 
lllOOl] . Still the uncertainty is large as compared to the extreme sensitivity of the fission lifetime. 
A prediction within four order of magnitudes is already great success 11931] . 


4.5. Correlation analysis for the key observables 

In a study of inter-correlations within the fission path (not shown here) we have found that 
the values of potential energy Vljif) along the fission path are highly correlated with each other 
which means, e.g., that knowing the barrier height By determines the fission potential V(J32) 
in a broad vicinity around the batTier. The same hold for the mass parameter. This allows to 
reduce the relevant information to the few key observables as indicated in figure [12] We will 
henceforth reduce the discussion to these. Eigure [15] shows the matrix of correlations between 
forces parameters and key observables of the fission path for SV-min. Remind that many of the 
SHE force parameters can be mapped to NMP which are a more intuitive way to characterise a 
force. The five bulk parameters of the SHE functional (volume term, density dependent term) can 
be represented by the five NMP equilibrium binding E/A with po, incompressibility K, symmetry 
energy asym and slope thereof dpasym- The two parameters for kinetic terms correspond to the two 
NMP for effective mass, m* / m and /ctrk- The two parameters for the gradient terms corresp ond t o 
the two surface energies Osurf and asurf.sym computed here in semi-classical approximation ill 12ll . 
Only pairing and spin-orbit parameters have no bulk equivalent and remain what they are. 

Amongst the key observables of the path (see figure [T2]l, we see a strong cori'elation between 
batTier By and curvature d^E(q)/dq^\^ , some cotTelation between these two quantities and the 
action integral, and very little cotTelation to the mass parameter at ground state. What the cor¬ 
relations of fission properties with force parameters is concerned, they are mostly small. Only 
m*/m and pairing parameters show some correlations, pairing particularly with the mass param¬ 
eter and m* Im with the barrier. After all, we see that there is no pronounced correlation of a 
fission observables to any particular force parameter. This indicates that each detail of the force 
contributes a bit to fission where from pairing is the strongest player. 


4.6. Trend analysis for the key observables o/^®®Hs 

Covariance analysis is a powerful tool to find out dependencies between observables (and 
force parameters). But it often lacks intuitive insight. As a complementing tool, one can perform 
a trend analysis. This is done by producing systematically various sets of fo rces w here in each 
set one forces parameter (or NMP) is varied. Such sets had been produced in starting with 
SV-bas as base point. The sets were obtained by constraining the same four NMP as for SV-bas 
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Figure 15: Matrix of correlations between the key observables of the fission path in ^®®Hs for selected force parameters 
SV-min. 


while varying one of the four NMP with respect to SV-bas. One can then draw the observables 
of interest versus the varied parameter. This makes the dependencies graphically visible which 
is often a helpful complementing illustration of the correlation. 

Figure[T6]illustrates the influences on fission in terms of the trends with the force parameters 
(NMP, respectively) K, m*fm, cisym, '<'trk, and pairing strength. In order to compare the trends 
within one panel (for each observable), we produced dimensionless force parameters. This is 
done by taking the parameters relative to the SV-bas point and scaled by the uncertainty from 
SV-min which characterises the allowed variation, i.e. the rescaled version of a parameter p 
is Prescai - (p - Psv-bas)/Apsv-min- This implies that the interval of rescaled force parameters 
[-1 , -Hi ] represents the range of reasonable parameters for which the r eproducti on of data remains 
nearly as good as for the minimum growing at most by one unit 11 lOl 111 , 921'). 

The trends in Figure [16] conhrm what we have seen already from the correlation matrix. 
Pairing strength has the strongest impact on fission. Next important is m* jm. There is also some 
impact from asym, but K as well as /ctrk play little role here in SHF. 


5. Summary and outlook 

We have presented a variety of theoretical results on hssion properties in super-heavy nuclei. 
Rather different modelling has been compared, the macroscopic-microscopic (MM) method, one 
the one hand, and self-consistent mean-held models on the other. The basis of the description 
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Figure 16: Trends of the key observables of the fission path of ^®®Hs with systematically varied force parameters. The 
parameters are drawn relative to their reference value from SV-bas and scaled to the extrapolation uncertainty of SV- 
min. These reference values and their uncertainties are given in the table above the plots. The values +1 for the scaled 
parameters stands for the variance shown in the table. The band ± 1 for the scaled parameters represents the allowed error 
band from analysis. 


is the same in both model families: the state of the system is characterised by a set of single¬ 
particle (SP) wavefunctions moving in a common mean-field and the hssion path consists in 
a series of such states with successively increasing deformation until two separated fragments 
emerge. The SP wavefunctions serve to account properly for the quantum shell effects which 
are crucial for binding of heavy elements and the topology of the potential energy surface. The 
fact that a couple of different mean-field states is involved means that the description deals, in 
fact, with collectively correlated states. These features are the same for all models. Different is 
the way in which the mean field and the total energy are dehned. The MM method models the 
mean-field directly as a (deformed) shell-model potential and composes the total energy from 
the macroscopic liquid-drop model plus shell corrections deduced from the s.p. energies. Self- 
consistent models start from an energy-density functional where-from the mean-field equations 
are derived variationally. The dedicated deformation which are needed to span the fission path 
are obtained by adding a (quadrupole) constraint. The density functionals manage to incorpo¬ 
rate a correct description of bulk properties (liquid-drop parameters) together with shell struc¬ 
ture. Within the self-consistent family, we have considered three different models, the Skyrme- 
Hartree-Fock (SHF) approach which is a non-relativistic density-functional method, the Gogny 
force which, again, uses non-relativistic kinematics but a finite-range force with a sprinkle of 
density dependence, and the relativistic mean-field model (RMF) which relies on a covariant 
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Lagrangian density functional. 

In spite of the very different modelling, the topological properties of the fission path are 
coming out the same way in all models: ground state deformations, fission isomers, shapes at 
barrier (reflection symmetry, triaxiality), islands of stability. These features are determined by 
quantum shell structure which is obviously correctly reproduced in all cases. Huge differences 
in the predictions are found in the quantitative aspects, height of fission barriers, and even more 
so for fission lifetimes. These values are very hard to control. Each detail of a model has an 
influence. Each one of the liquid-drop bulk properties has some impact on barriers, the effec¬ 
tive mass which regulates shell effects plays also a role, and last not least pairing is a crucial 
ingredient. Statistical analysis on top of least-squares fits within the SHE approach estimates the 
uncertainties on barriers typically of the order of +1 MeV. We have to add substantial leeway 
for systematic errors which, however, cannot be estimated in systematic manner. The situation 
is even more involved for fission lifetimes. Eurther ingredients are required for its computation: 
entry point for spontaneous fission, collective mass along the path, and quantum corrections for 
spurious zero-point energies. All these are extremely subtle quantities which depend very sensi¬ 
tively on all details of a model and on the way they are computed. There is lot of development 
work yet to be done to reach better agreement in the quantitative predictions. Eor the time being, 
the method of choice is to calibrate barriers and lifetimes for the heaviest known nuclei to allow 
safer predictions for the even heavier elements. 

Most self-consistent calculations describe fission as a one-dimensional, adiabatic process. 

MM calculations indicate that this may not always be sufficient. There are situations in which 
deformation valley and fission valley go in parallel for while. It is likely that one needs here a 
description in terms of multi-dimensional collective motion. This is an extremely challenging 
problem which has not been attacked thoroughly yet. We hope that increasing computing power 
allows steps into that direction. The quality of the adiabatic approximation deserves also critical 
testing. There exist already investigations on dynamical aspects of fission. This line of develop¬ 
ment needs also further inspection to learn more about the limits of presently used methods and 
possibly to improve on them. 

Self-consistent models are optimized by least-squares fits to data. Trend- and correlations- 
analysis built on the underlying least-squares fits shows that the relations of fission properties 
and forces parameters and NMP are mixed. We see the reduction of uncertainties when fixing a 
couple of NMP. But there is no single NMP which has predominant influence. Eission properties 
depend on every detail of the parametrisation, the bulk properties given by NMP as well as 
pairing strength and spin-orbit coupling (determining shell effects). 

All this shows that a proper description of fission is extremely demanding. More research is 
needed to achieve a higher reliability for predictions of fission lifetimes in the context of self- 
consistent mean-field models. 
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